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Abstract 

Cosmological simulations of disc galaxy formation, when compared to the observed Tully-Fisher relation, 
suggest a low Mass-to-Light (M/L) ratio for the stellar component in spirals. We show that a number of 
"bottom-light" Initial Mass Functions (IMFs) suggested independently in literature, do imply M/L ratios as 
low as required, at least for late type spirals (Sbc/Sc). However the typical M/L ratio, and correspondingly 
the zero-point of the Tully-Fisher relation, is expected to vary considerably with Hubble type. 

Bottom-light IMFs tend to have a metal production in excess of what is tipically estimated for spiral 
galaxies. Suitable tuning of the IMF slope and mass limits, post-supernova fallback of metals onto black 
holes or metal outflows must then be invoked, to reproduce the observed chemical properties of disc galaxies. 
Keywords: Galaxies: spirals, chemical and photometric evolution; stars: Initial Mass Funtion 
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Figure 1: Straight lines: observed Tully-Fisher rela- 
tion for Sbc-Sc disc galaxies (Dale et al. 1999; h=0.7 
adopted here), for different assumptions about the 
stellar M/L/ ratio. Triangles: simulated galaxies by 
Sommer-Larsen et al. (2003). Also shown is the loca- 
tion of the Milky Way and NGC 2841. 

1 Introduction 

Recent N-body+SPH cosmological simulations of the for- 
mation of disc galaxies reproduce the observed Tully-Fisher 
(TF) relation (Dale et al. 1999), provided the mass-to- 
light (M/L) ratio of the stellar component is rather low, 
M/L/=0.7-l in the I-band (Sommer-Larsen & Dolgov 2001; 
Sommer-Larsen et al. 2003; Fig. 1). The location of the 
simulated galaxies in the (M», V c ) plane of Fig. 1 is quite 
independent of the adopted Initial Mass Function (IMF) 
or feedback efficiency: the baryonic mass that cools out 
to form a galactic disc and its resulting circular velocity 
correlate so that data points tend to move along the TF 
relation, hardly affecting the zero-point (Navarro & Stein- 
metz 2000ab). However, the IMF is crucial for the M/L ra- 
tio, to translate the stellar masses M, to luminosities and 
compare the simulated TF relation to the empirical one. 



Although the zero-point of the simulated TF may change 
with the concentration of the dark matter halos, and hence 
with the normalization of the power spectrum as (Navarro 
& Steinmetz 2000ab; Eke et al. 2001), many other argu- 
ments support a low stellar M/L ratio in spiral galaxies: 

The stellar mass of the Milky Way is M* ~ 5 x 10 10 M ; 
to lie on the observed TF relation as other spirals, its M/L/ 
must be < 1 (Sommer-Larsen & Dolgov 2001; Fig. 1). A 
low M/Lj < 0.8 is also derived for the massive Sb galaxy 
NGC 2841, when compared to the observed TF relation 
(Portinari et al. 2004a, hereinafter PST; Fig. 1). 

Based on bar instability arguments, Efstathiou et al. 
(1982) suggest an upper limit of M/L_g < 1.5 h for discs, 
i.e. M/Lb < 1 for h=0.7. (h indicates the Hubble constant 
Ho in units of 100 km sec" 1 Mpc~ ). 

The stellar M/L ratio is related to the issue as to whether 
discs are maximal or sub-maximal, i.e. as to whether they 
dominate or not the dynamics and rotation curves in the 
inner galactic regions. Even in the case of maximal stellar 
discs, lower M/L ratios for the stellar component are re- 
quired, than those predicted by the Salpeter IMF (Bell & 
de Jong 2001). And it is still much debated whether discs 
are maximal or sub-maximal; for his favoured sub-maximal 
disc model, Bottema (2002) finds M/L/ ~ 0.82. 

Finally, two recent dynamical studies of individual spi- 
ral galaxies yield M/L~l in the B, V and I band for the Sc 
galaxy NGC 4414 (Vallejo et al. 2002) and M/L/=l.l for 
the disc of the Sab spiral 2237+0305, Huchra's lens (Trott 
& Webster 2002). 

In this paper we discuss if M/L ratios so low are com- 
patible with our understanding of stellar populations and 
chemical evolution in disc galaxies. We also address the ef- 
fects of different star formation histories on the TF relation 
for different Hubble types. 



2 Star Formation History 
and Initial Mass Function 

The M/L ratio of the stellar component of a galaxy (includ- 
ing both living stars and remnants) depends on the stellar 
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Figure 2: Mean 6-parameter as a function of Hubble 
type. Data from Sommer-Larsen et al. (2003). 

Initial Mass Function (IMF) and on the star formation his- 
tory (SFH) of the system. 

The SFH of a disc galaxy is related to its Hubble type: 
Kennicutt et al. (1994) demonstrated that the sequence of 
spiral types is in fact a sequence of different SFHs in the 
discs, as traced by the birthrate parameter 

SFR 

~ < SFR > 

or the ratio between the present and the past average star 
formation rate (SFR); see Fig. 2. Our reference TF relation 
by Dale et al. (1999) in Fig. 1 is representative for Sbc-Sc 
spirals (Giovanelli et al. 1997). Kennicutt et al. find that 
Sbc-Sc discs correspond to b=0.8-l, hence we will concen- 
trate our discussion on objects with SFHs corresponding to 
these values for 6. 

As to the IMF, a standard Salpeter slope extended 
over the mass range [0.1-100] Mq certainly yields much 
higher M/L ratios than those mentioned in the introduction 
(Fig. 4, top left). There is however plenty of observational 
evidence that the IMF flattens below ~ 1 Mq , possibly with 
a turn-over at low masses, and hence is "bottom-light" 
with respect to the Salpeter IMF. 

We considered the following IMFs: the Salpeter (1955) 
IMF; the Kroupa (1998) IMF, derived from field stars in 
the Solar Neighbourhood; the Kennicutt IMF, derived 
from the global properties of spiral galaxies (Kennicutt 
et al. 1994); the Larson (1998) IMF, with an exponential 
cut-off at low masses as favoured by recent determinations 
of the local IMF down to the sub-dwarf regime (Chabrier 
2001, 2002); a modified Larson IMF, with the same peak 
mass but a steeper slope at high masses, in accordance with 
local determinations (Scalo 1986); the Chabrier (2001, 
2002) IMF, derived from low mass stars and brown dwarfs 
in the local field. With respect to Salpeter, all the other 
IMFs are "bottom-light" (Fig. 3). 



3 Results from simple models 

We will focus on the M/L ratio in the I-band, which is an 
excellent optical tracer of stellar mass, because the I-band 
luminosity is less sensitive to recent sporadic star formation 
and to corrections for dust extinction than bluer bands. 
Besides, the M/L ratio of Single Stellar Populations (SSPs) 
in the I-band is less metallicity dependent than in other 



Figure 3: Comparison between the Salpeter IMF and 
the "bottom-light" IMFs considered in this paper, all 
normalized to the same total integrated mass of 1 Mq . 

bands, and it is less sensitive to the specific treatment of 
the AGB phase than redder, NIR bands (PST). 

We computed SSPs based on the latest Padua isochrones 
(Girardi et al. 2002), for the different IMFs in §2. These are 
to be convolved with suitable SFHs. We generated a set of 
SFHs characterized by different values of the fe-parameter, 
by adopting an exponentially decaying SFR oc e~~ with 
different decaying rates. SSP metallicities around solar can 
be considered typical for spiral galaxies. The resulting M/L 
ratio of the global stellar population (including remnants) 
as a function of the ^-parameter is displayed in Fig. 4, for 
the different IMFs. The range 0.8 < b < 1 is representative 
of Sbc-Sc spirals. The red shaded area indicates the range 
in M/L=0.7-l suggested in §1. 

Fig. 4 shows that bottom-light IMFs can in fact yield 
M/L/ < 1 for late-type spirals, though one probably needs 
slightly "lighter" IMFs than the local Kroupa one. 



4 Offsets of the TF relation with Hubble type 

From Fig. 4, the stellar M/L ratio is expected to vary with 
Hubble type due to the differences in SFH parameterized by 
b. This effect implies systematic offsets with Hubble type of 
the luminosity zero-point of the TF relation (Rubin et al. 
1985; Giovanelli et al. 1997; Kannappan et al. 2002). 

Fig. 5 shows the M/L ratio as a function of b, normal- 
ized to the value corresponding to 6=1. The scale on the 
right axis indicates the corresponding shift in magnitude. 
With respect to Sbc/Sc spirals, we predict a systematic off- 
set of 0.3-0.4 mag for Sb's (b ~0.35) and of 0.6-0.8 mag for 
Sa/Sab's (b ~0.1), as a result of the different characteristic 
SFHs. These offsets are only slightly reduced when bulges 
are added to discs in the computation of the total M/L 
ratios of galaxies (PST). The offsets we predict are larger 
than the empirical ones found by Giovanelli et al. (1997): 
0.1 mag for Sb spirals and 0.32 mag for earlier types. How- 
ever, the extent of the observed offsets in the TF relation 
is still a matter of debate: for instance, the larger offsets 
found by Kannappan et al. (2002) in the R-band, 0.76 mag 
for Sa's, are in good agreement with our predictions. 

The characteristic SFH of a galaxy is traced by its 
colours, so that the offsets in M/L ratio due to different 
SFHs can be corrected for, by applying suitable M/L vs. 
colour relations (Bell & de Jong 2001; PST). 



M/L ratio, IMF and chemical evolution in disc galaxies 
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Figure 4: I-band M/L ratio at varying ^-parameter of the SFH, for different IMFs. The red shaded area marks 
the range M/L/=0.7-l favoured by observations for Sbc-Sc discs (corresponding to 6=0.8-1). 
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Figure 5: Relative M/L ratio normalized to the value 
corresponding to 6=1 models. 

5 Bottom light IMFs and chemical evolution 

Besides simple models with exponentially declining SFHs, 



confirmed: the "bottom-light" IMFs considered here im- 
ply M/L/ <1 for Sbc/Sc spirals, as required in §1. This 
conclusion remains valid also when the contribution of the 
bulge to the global M/L ratio is included (PST). 

Chemo-photometric models allow also an insight on 
the implications of "bottom-light" IMFs for chemical evo- 
lution. Some of the IMFs considered (Kennicutt, Larson 
and Chabrier) are too efficient in metal production to re- 
produce the observed properties of spirals. In particular, 
the resulting gas fraction is much larger than the observed 
Mg a3 /LB ~ 0.5 Mq/Lq (solid dots in Fig. 6, bottom panel). 
For a given IMF and corresponding metal production, the 
final gas metallicity predicted by chemical models increases 
at decreasing gas fraction (Tinsley 1980, Pagel 1997). Our 
models, calibrated to reproduce the observed metallicities, 
tend to predict too high gas fractions. Conversely, if they 
were calibrated to reproduce the observed gas fractions, 
they would result in too high metallicities. 

This excessive metal production is readily understood 
since the enrichment efficiency of a stellar population, or 
its "net yield", is inversely proportional to the mass frac- 



we also computed more realistic, multi-zone chemo-photometric 

tion that remains forever locked m low-mass stars and rem 



models of galactic discs, including infall, inside-out forma- 
tion and radially varying star formation efficiency. Chem- 
ical evolution is followed with the code by Portinari et al. 
(1998), Portinari & Chiosi (1999) and the models are cali- 
brated to reproduce the typical metallicity and metallicity 
gradient of Sbc/Sc discs (PST) . Six sets of models have been 
computed for the six IMFs in §2. The corresponding pho- 
tometric properties are calculated by convolving the SFH 
and metal enrichment history of each annulus of the disc, 
with a grid of SSPs of metallicities between 5 x 10~ 4 Z@ and 
5Z Q (PST). 

For each IMF, the chemo-photometric models confirm 
the M/L ratios predicted, as a function of b, by the simple 
models in Fig. 4; see the example for the Kennicutt IMF 
models in Fig. 6, top panel. The results of §3 are thus 



nants (Tinsley 1980; Pagel 1997); and for "bottom-light" 
IMFs the locked-up fraction is small. This effect can be 
compensated by a steep slope at the high-mass end, which 
reduces the number of massive stars and the related metal 
production: the Kroupa or modified-Larson IMFs, for in- 
stance, with a steep Scalo slope do not overproduce metals 
(see PST for details) . A steep slope for the integrated field 
stars IMF is expected, in fact, if stars form in star clus- 
ters of finite size even from an intrinsically shallower IMF 
(Kroupa & Weidner 2003). 

For those bottom-light IMFs that do imply an excess in 
metal production, one way to reconcile models with obser- 
vations is to tune the upper mass limit of the IMF, again 
reducing the number of massive stars and related metal 
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Figure 6: Top panel: M/L ratio for chemo-photometric 
models with the Kennicutt IMF (solid dots) ; they con- 
firm the prediction from the simple exponential mod- 
els (black line) and are within the observed range of 
M/L/=0.7-l (red shaded area). Bottom panel: Pre- 
dicted gas fractions compared to the observed range, 
marked by the red dashed lines. Green triangles: mod- 
els with IMF mass limits tuned between 0.05-0.1 M Q 
and 30-35 M Q , rather than the standard mass range 
[0.1-100] Mq, to reproduce the observed gas fraction. 

production — or equivalently assume that above a certain 
progenitor mass, the metals produced by a supernova fall 
back onto a black hole after the explosion. The observed gas 
fractions can then be matched without altering the stellar 
M/L ratio significantly (green triangles in Fig. 6). 

Alternatively, we need to invoke substantial outflows of 
metals from disc galaxies into the intergalactic medium, to 
reconcile the high enrichment efficiency with the observed 
low gas fractions. This behaviour is reminiscent of that of 
elliptical galaxies, responsible for the enrichment of the hot 
gas in clusters. 

With "standard" IMFs suited to model the chemical 
evolution of the Solar Neighbourhood (e.g. the Kroupa IMF) 
it is impossible to account for the observed metal enrich- 
ment in clusters (Portinari et al. 2004b). Possibly, some 
of the "bottom-light" IMFs advocated here to reproduce 
low disc M/L ratios, suggest a scenario where the IMF 
and the enrichment efficiency may be the same in spiral 
and cluster galaxies, and in both cases much of the metals 
are dispersed into the intergalactic medium. However, sub- 
stantial outflows would challenge our understanding of disc 
galaxy formation: in galactic discs, star formation proceeds 
at a smooth, non burst-like pace and the observed "foun- 
tains" and "chimneys" do not have enough energy to escape 
the galactic potential; winds are far less plausible than in 
spheroids. Moreover, strong ongoing stellar feedback and 
outflows could significantly hamper the dynamical forma- 
tion of galactic discs from the cool-out of halo gas. 

The alternative to major winds from disc galaxies is 



that the metal production in spirals is much lower than 
in galaxy clusters, because of a different IMF (Portinari 
et al. 2004b). The IMF may change out of Jeans-mass 
variations with redshift (Moretti et al. 2003 and references 
therein); or be a universal function within star clusters, 
but lead statistically to more high-mass stars in massive 
ellipticals, where in regimes of intense star formation larger 
star clusters can be formed (Kroupa & Weidner 2003). 
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